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1. Introduction 



1.1. Equidistribution theorems for closed geodesies. Let X be a com- 
pact surface with a metric of constant negative curvature k = —1, SX be 
the unit tangent bundle of X, and : SX SX the geodesic flow. We 
think of SX as the set of initial conditions {z, Q with z & X the position 
and ( the direction vector. 

The geodesic flow is ergodic with respect to Liouville measure dx, the 
smooth invariant measure for the flow: Generic geodesies become cquidis- 
tributed, in the sense that for Lebesgue-almost all initial conditions xq G 



for integrable observables on SX. 

As is well known , there are inflnitely many closed geodesies, in fact the 
number 7r(T) of closed geodesies of length at most T grows exponentially 
with T: 7r(r) ~ e'^/T as T ^ oo [40], [8]. For a closed geodesic C, let i{C) 
be its length and nc be the arc-length measure along C, i.e. for F G C{X), 



This is a measure on SX, invariant under the geodesic flow and of total 
mass £{C). Closed geodesies become, on average, uniformly distributed 
with respect to dx-} For any observable F G C{SX) we have 



Lallcy [20] determined the fluctuations of the numbers fic{F)/ y^£{C) for 
F as above of zero mean, as C varies over closed geodesies ordered by length. 
He showed that they are Gaussian with mean zero and variance V{F,F) 
where V is the hermitian bi-linear form on functions of zero mean given by 



In variable negative curvature, one needs the Bowen-Margulis measure here. To get an 
equidistribution statement involving Liouville measure, one needs to weigh each geodesic 
by its "monodromy" . 



SX, 






(1.1) 
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The negative curvature guarantees that the correlations in the inner integral 
decay exponentially as t ^ ±00, so that V is convergent [33]. 

The bilinear form V is positive semi-definite, and is degenerate, in fact if 
Fq is smooth then V{Fq, F) = for all -F if and only if Fq is a derivative in 
the flow direction: Fq = ^|i=o/o$* for some other observable / e C'^{SX). 

An important involution of SX is time reversal symmetry 

w : {z,C) ^ {z,-0 

which reverses the direction vector of the initial condition, and satisfies 
w o $* = $~* o w. It induces an involution on the set of geodesies, taking a 
geodesic C = {$*a;o : t G M} to its time reversal C = -wC = {<I>*wxo : s G 
M}. 

Time reversal symmetry can also be incorporated in Lalley's theorem: 
To do so, note that for a closed geodesic C, its time-reversed partner C is 
also closed and both have the same length: £{C) = £{C). Grouping these 
together yields the measure /x^^" := /ic + mass 2i{C). By Lalley's 

theorem, the fluctuations of ^'^^^ / \J2(.[C) arc again Gaussian with mean 
zero but with variance given by the hermitian form 

(1.2) y^^'='^(Fi, F2) = V {Fl"^'', Fl"^"") 

where F"^^^"- = {F + F o ■w)/2 is the even part of F under w. Note that 
^gjen jg invariant and V^^^"' is bi-invariant under the geodesic flow as well 
as under time-reversal symmetry w. Both hermitian forms V and V^'"^'^ on 

LliSX) := {/ G L\SX) : / fix)dx = 0} 

Jsx 

can be diagonalized and computed explicitly by decomposing the regular 
representation of PSL2(M) on this space, see § 3. 

1.2. The modular surface. In this paper we investigate fluctuations of 
measures on the modular surface associated with grouping together geodesies 
of equal discriminant. As is well known, any of our compact surfaces X 
may be uniformized as a quotient of the upper half-plane H, equipped with 
the hyperbolic metric, by a Fuchsian group F. Furthermore, the group 
G = PSL2(M) of orientation preserving isometries of H acts transitively on 
the unit tangent bundle SX, giving an identification SX ~ F\G, reviewed in 
§ 2. The modular surface is obtained by taking F = PSL2(Z); the resulting 
surface is non-compact (but of finite volume) and has elliptic fixed points, 
but these issues will not be important for us. 

Closed geodesies correspond to (hyperbolic) conjugacy classes in F, with 
the length of a closed geodesic C given in terms of the trace t of the cor- 
responding conjugacy class by i{C) = 21og(t + VW^)/2. In the case of 
the modular surface, the hyperbolic conjugacy classes correspond to (strict) 
equivalence classes of integer binary quadratic forms ax^ + bxy + cy^ (also 
denoted by [o, 6, c]), of positive discriminant d := — Aac with the modu- 
lar group acting by linear substitutions (we need to exclude discriminants 
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which are perfect squares). The discriminant disc(C) of a closed geodesic C 
is defined as the discriminant of the corresponding binary quadratic form. 

For d > 0, d = 0, 1 mod 4 and d not a perfect square, let /i, . . . fuid) be 
the classes of binary quadratic forms of discriminant d. We do not assume 
that fj = [aj.,hj^Cj\ is primitive and so H{d) is the Hurwitz class number 
[21]. Let 

(1-3) ed = ^ , id>0,tid>0 

be the fundamental sohition of the Pellian equation — du? = 4. Then as 
in [37, 39] associate to each fj the F-conjugacy class (it is well-defined) of 
the matrix 



(1.4) 



2 



This gives H(d) closed geodesies for each discriminant d, all of length 2 log e^. 
Let be the corresponding measure on SX: 

(1-5) = XI ■ 

disc{C)=d 

These measures are the arithmetic measures in the title of the paper. 
They have been studied extensively and the primary result about them is 
that they become equidistributcd as d ^ oo. That is, if F is bounded and 
continuous on SX and has mean zero, then 

0, as a — DO . 



H{d)2\oged 

Linnik [24] developed an ergodic theoretic approach to this equidistribu- 
tion problem and recently [3] have shown that this method leads to a proof 
of this specific result. The first proof of equidistribution is due to Iwaniec 
[11] and Duke [1]. Iwaniec established the requisite estimate for Fourier 
coefficients of holomorphic half-integral weight forms (of weight > 5/2) and 
Duke obtained the estimates for weight 3/2 and weight zero Maass forms. 
In view of our reductions in Sections § 3 and § 6, these together imply^ the 
full equidistribution on SX. 

The measures jid enjoy some symmetries (see [39]). Firstly they are in- 
variant under time reversal symmetry: w^u^ = /i^. Secondly, let r be the 
involution of F\G given by g i— > 5~'^g5, where 



5 = 



1 
-1 



Specifically, by (3.25), (6.1) and (6.3), equidistribution on SX is reduced to estimation 
of Fourier coefficients of classical holomorphic forms of half integer weight and Maass forms 
of weight 1/2. 
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(it is well defined since 6 ^Td = T). In terms of the coordinates (zX) on 
SX, r is the orientation-reversing symmetry 

r : {z, C) ^ i-z, -C) ■ 

The measure /x^ is also invariant under r. The involutions w, r commute 

and their product rw is also an involution. Thus fi^ is invariant under 
the Klein four-group H = {/, r, w, rw}. These involutions induce linear 
actions on L'^{T\G) by /(x) i-^ fiHx)) with h e H and x G r\G and 
we denote these transformations by the same symbols. The fluctuations 
of the measures fi^ inherit these symmetries and since we are particularly 
interested in comparing their variance with the classical variance V we define 
the symmetrized classical variance V^'^"^ on functions of mean zero on r\G 

by 

(1.6) V'y"'{Fi,F2) := ViFl^"", "*) 



where 



psym _ 1 _ 

4 ^ 

heH 



1.3. Results. We can now state our main results about the fluctuations of 
lid- We normalize these measures as 

This is essentially equivalent to normalizing by the square root of the total 
mass, ■sjH(d)2\og e^, see Remark 1.4.2. The space of natural observables 
for which one might compute these quantities is LQ(r\G), or at least a dense 
subspace thereof. This space decomposes as an orthogonal direct sum of the 
cuspidal subspace 

^Lp(r\G) := {/ e L\T\G) : / f{nx)dn, for a.e. x e r\G} 

JNr\r\N 

where N = {(^ : u G M}, and the unitary Eisenstein series [5]. The 

former is the major and difficult part of the space LQ(r\G) so we will con- 
centrate exclusively on it. One can easily extend our analysis of the variance 
to the unitary Eisenstein scries. 

Theorem 1.1. Fix smooth, K-finite Fi,F2 G {T\G). Then 



(1.7) lim 



-'cuspy 



and there is a limiting variance 



(1.8) BiF,,F2)=li^^-—-—-Y: 



Y^oc #{d:d<Y} ^ dV4 ciV4 
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We call this variance B the "arithmetic variance" . The structure of the 
bilinear form B is revealed by choosing a special basis of observables, com- 
patible with the symmetries of the problem. Recall that the unit tangent 
bundle is a homogeneous space for G = PSL2(M), and thus it is natural to 
decompose the space L^(r\G) into the irreducible components under the 
G-action. In addition, there is an algebra of Hecke operators acting on this 
space, commuting with the G-action, hence also acting on each isotypic G- 
component. We take observables lying in irreducible spaces for the joint 
action of G and the Hecke operators - the automorphic subrepresentations 
of L^„^p(r\G). Denote the decomposition of the regular representation on 
Ll^gp{T\G) into G- and Hecke-irreducible subspaces by 

oo 

(1-9) Ll,spiAG) = ^W^. 

3=1 

SO TTj is a cuspidal automorphic representation. 

In order to describe the arithmetic variance explicitly we need a more 
detailed description of the W-j^.^s. To each ttj is associated an even integer k, 
its weight (see § 3) which we indicate by tTj. For k = there are infinitely 
many tt^'s corresponding to Hecke-Maass cusp forms on X, while for A; > 
there are such tTj (where is cither [^] or [^] + 1, depending if k/2 = 1 
mod 6 or not), corresponding to holomorphic Hecke cusp forms of weight k. 
For A; < let 

W^k = = {J : f e W-k} 

for j = 1,2, and these correspond to the anti- holomorphic Hecke 

cusp forms. With these we have the orthogonal decompositions 

oo rffc 

j=l k>12 j=l ^ 

oo (ifc 

(1-10) = E^-°®EEf^^^ 

j=l k>l2 j=l 

where 

(1.11) U^o = W^o, and U^k = W^k © W^-k 

To each ttj as above one associates an L-function L{s,-Kj) given by 

oo \ / \ 

(1.12) L(,,^^.) = ^^^, 3ft(,)>l 

n=l 

where (n) is the eigenvalue of the Hcckc operator T„ acting on W^.. It 
is well known (Hecke-Maass) that L{s,TTj) extends to an entire function and 
satisfies a functional equation relating its value at s to 1 — s. In particular 
the arithmetical central value L{^,TTj) is well defined (and real). 
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Theorem 1.2. Both and B are diagonalized by the decomposition 

(1.10) and on each subspace U k we have that 

j 

i ^ i 

where c(0) = G/tt and c{k) = I/tt if k > 0. 
1.4. Remarks. 

1.4.1. The hermitian forms V^'^^ and B can be computed exphcitly on each 
U^k (sec § 3). Time-reversal symmetry w forces to vanish on U^k for 
k = 2 mod 4. Also orientation-reversal symmetry r fixes the weight zero 
spaces U a and hence takes the generating vector (see ^ 3) (j)^ €z ir^ into 

^(pj- Corresponding to this sign we call U^o even or odd. According to § 3, 

ysym jg completely determined on U k by its value on the generating vector; 

j 

hence it follows that V^y^lu^ = for the odd 7r?'s. In the above cases where 

I TVj J 

y^^"^\u fc vanishes, the sign e-j^j of the functional equation of L{s,Trj) is — 1 

and hence the central L- value L(^, tTj) = for reasons of symmetry. In the 
other cases {k = mod 4 and 7r° even), ejr^ = 1 and V^y"^\u„^ 7^ 0. One 
expects that in these cases L{^,'Kj) 7^ as well. However if we pass from 
r = PSL2(Z) to a congruence subgroup, where our analysis can be carried 
over with similar results, then there will be tt's corresponding to holomorphic 
forms for which L(^,7r) = for number-theoretic reasons, specifically the 
conjecture of Birch and Swinnerton-Dycr [43] . In this case the restriction of 
the arithmetic variance to such a subspace will vanish for reasons far deeper 
than just symmetry. 

1.4.2. The normalization iX(i{F)/d^f^ is natural from the arithmetic point of 
view. To be consistent with the previous normalization we should use the 
square root of the total mass ^yH{d)2 log e^ of the measure. By Dirichlet's 
class number formula for d fundamental, when H(d) = hid) is the ordinary 
class number (and similar formulae for all d), 

(1.13) h{d)\oged = ^fdL{l,Xd) . 

The fluctuations of L(l, Xd) are mild and well-understood [4] and hence the 
normalizations are essentially the same. In any case one could use methods 
as in [12, Chapter 26] to remove the weights L(l, Xd) and deduce Theorem 1.1 
with this other normalization. 

1.4.3. In § 3 we show, in a more abstract context, that the space of lin- 
ear forms on an irreducible unitary representation of G which are invari- 
ant by both the geodesic flow and time reversal symmetry is at most one- 
dimensional, and how to incorporate orientation-reversal symmetry. This 
shows that the form that the arithmetic and "classical" variance take is uni- 
versal. That is for any family of such invariant measures, the variance B\ 
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if it exists, is determined completely in each irreducible representation of G 
by B'iyo^vo), where vq is either a spherical vector, or a lowest (or highest) 
weight vector in the representation. 

1.4.4. The geometric problem is to order the by the length of any of 
the geodesic components of the measure. We do not know how to do this. 
What we can do is to compute the variance of the //^'s when ordered by the 
discriminant d. Prom the arithmetic point of view this ordering is anyway 
the most natural one. For many considerations these two orderings of 
yield quite different answers (see [39]). However for the fluctuations we 
believe they are similar. 

The difficulty in proving the same result of the /x^'s ordered by td (or e^) 
is apparent already for Fi = F2 = f a holomorphic cusp form of weight 
m = mod 4. In this case according to the formula of Kohnen and Zagier 
[17], we have for d a fundamental discriminant say 

(1.14) ^ = *L(^,/®x.) 

(with * explicit and under control). Thus we would need to understand the 
averages 

(1-15) ^L(i/®xd) 

td<Y 

The first, but big, step in this direction would be to understand 
(1-16) 

t<Y ^ 

(see [34] for an execution of such an analysis on a simpler problem). This 
appears to be beyond the well developed techniques for averaging special 
values of L-functions in families. We leave it as an interesting open problem. 

1.4.5. The recent work [36] giving lower bounds for moments of special values 
of L-functions in famihes, together with (1.14), shows that the fluctuations 
of ^d[F)/d^^^ are not Gaussian, at least not in the sense of convergence of 
moments. 

1.4.6. The arithmetic variance B in Theorem 1.1 is the same as the quan- 
tum variance for the fluctuations of high energy eigenstates on the modular 
surface that were calculated in [25] and [44]. We expect that the variance 
for the iidS when ordered by length will be the same as B. This would give 
a semi- classical periodic orbit explanation for the singular finding [25] that 
the quantum variance is B rather than V^'"^^. It points yet again, just as 
for the local spectral statistics (see the survey [38]), to the source of the 
singular behaviour of the quantum fluctuations in arithmetic surfaces being 
the high multiplicity of the length spectrum. Similar phenomena are found 
for the quantized cat map [18, 19]. 



VARIANCE OF ARITHMETIC MEASURES ON THE MODULAR SURFACE 9 

1.5. Plan of the paper. We end with an outline of the paper and the 
proof of Theorem 1.1. In § 2 we give recall some background connecting the 
dynamics on the modular surface with the group structure on SL2(IR). In 
§ 3 we show that up to a scalar multiple, there is at most one linear form on 
the smooth vectors of an irreducible unitary representation of SL2 (M) which 
is invariant under the action of the diagonal subgroup (corresponding to 

the geodesic flow) and the element ^ corresponding to time-reversal 

symmetry. We show that such a linear form is determined by its value on 

a "minimal" vector - a spherical vector in the case of a principal series 
representation and a lowest /highest weight vector for holomorphic/anti- 
holomorphic discrete series representations. We then bring in invariance 
under orientation reversal and apply the results to show that the bilinear 
forms V^y^ and B are determined by their values on Maass forms and holo- 
morphic modular forms. 

In §4 we give present some background on half-integral weight forms, and 
in § 5 we discuss Rankin-Selberg theory for these, giving a mean-square 
result for Fourier coefficients along positive integers by modifying work of 
Matthes [28] for weight zero forms. 

In § 6 wc review the results of Maass [26], Shintani [41], Kohnen [15, 16] 
and Katok-Sarnak [14], relating periods along closed geodesies to Fourier 
coefficients of theta-lifts. This allows us to express //^(-F) in terms of Fourier 
coefficients of half-integral weight forms on ro(4); the precise normalizations 
in terms of the inner products of the forms and their 0-lifts arc crucial here. 
This is where the factor L(i,7r) appears. These results put us in a position 
to use the Rankin-Selberg theory of § 5 to determine the variance which 
we do in § 7. 

1.6. Acknowledgments. We would like thank Akshay Venkatesh for in- 
sightful discussions of the material related to this paper. Supported by 
NSF FRG Grant DMS-0554373 (Sarnak and Luo) and by the Grant No 
2006254 (Sarnak and Rudnick) from the United States-Israel Binational Sci- 
ence Foundation (BSF), Jerusalem, Israel. 

2. Background on periods 

2.1. The upper half-plane and its unit tangent bundle. We recall the 
hyperbolic metric on the tangent bundle of the upper half-plane H = {z = 
X + iy : y > Q}. We identify the tangent space at 2; G EI with the complex 
numbers: T^H 2± C. The hyperbolic metric on T^H is then given by 

and the unit tangent bundle S'H is then identified with 



{(z,OeHxC:|C| = $>(z)} 
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2.1.1. Isometries. A unimodular matrix g = {'^ j J ^ SL2(M) acts on the 



c d J 

upper half-plane H via z ^ {az + b)/ {cz + d) . Set 

3{g,z) =cz + d 

The differential of the map is g'{z) = {ad — bc)/{cz + d)"^ = {cz + = 
1/ j{g,zY . The induced map on the tangent bundle TH is then 

{z,i)^{9{z),g'{z)0 

Note that this is an action: \i g^h ^ SL2(M) then g{h{z, Q) = {gh){z, Q. A 
computation shows that we get an isometry of H: 

{^,v)z = {g'{z)^,9'iz)ri)g(^^) 

2.1.2. Group theory. Define matrices 

n(x)-(^ ^^ a(y)-(y'^' \ _ r cos(</</2) sin(0/2)\ 

""^^^ - V V ' " V y-'/' J ' ^ ' ~ [-MH'^) cos((^/2) J 

The rotation K{(f)) preserves the base point i = G H. Note that 

k{(I) + 2tt) = -K{(f>) 

and thus we get the same element in PSL2(M). 
Setting g^^y^^ = n{x)a{y)K{(f)) we find 

9x,yA^^^) = {x + iy,iye^^) 

so that using the basepoint {i, i) G SM. of the upward pointing unit vector 
at i = \f^-i G H, we get a bijection 

PSL2(M)~SM, g^g{i,i) 

We may then identify functions on PSL2(M) and on SM: If F{z,Q is a 
function on SM we may define F on SL2(M) by 

F{g) := F{g{i,i)) 

so that F{gx^y^^) = F{x + iy, iye"t'). 

2.1.3. Geodesies. The geodesic flow on SM. is deflned by : {z,Q i— >■ 
{z{t),({t)) being the endpoint of the (unit speed) geodesic starting at z 
in direction C = iye^'^. It turns out that on PSL2(M) the geodesic flow is 

multiplication on the right by I _^/2 , that is 



e 



-t/2 

) ih^) 



Indeed, for an initial position {z, C.) G 51HI, we write {z, C) = g{i, i) and the 
the geodesic 7(t) = $*(2;, C) starting at (^;, C) will be the translate by gx,y,^ 
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of the geodesic 70 (i) starting at the initial condition {i, i): ^(t) = gx,y,<f>"yo{t)- 
A computation shows that 



70 (t) = ie'i,eH) = ( ^ _j/2 ) 

and therefore 



lit) = gx,y,.l> { g_t/2l(^,*) 



2.1.4. Time-reversal symmetry. A fundamental symmetry of phase space 
SM is time reversal {z,Q '-^ i^^—O- Using it, one has a symmetry of the 
set of geodesies, corresponding to reversing the orientation. In PSL2(M) it 

is given as g >—>■ gw, where = ■ Indeed, if {z, Q = g{i, i) G SM 

then 



{z,-0 = g{i,-i) = g (^^ 



2.1.5. Orientation reversal. Another fundamental symmetry is orientation 
reversal {z, () 1— (^—z, —(). On PSL2(M) it is given by the map 



g 6g6, S 



1 
-1 



2.1.6. K -types. Let k be an integer. Suppose that F : SM — C satisfies 

Fiz,e"'0 = e'^''F{z,0 
Then the corresponding function F on PSL2(M) satisfies 

FigK{a)) = e'^''F{g) 

that is transforms under under the right action of the maximal compact 

K = S'0(2)/{ib/} via the character K{a) ^ e*'^". As an example we start 
with a function / on H and define Ff{z, Q = C^f{z). 

2.2. Quotients. Let T C PSL2(M) be a Fuchsian group , M = r\H and SM 
the unit tangent bundle to M. The identification SM. ~ PSL2(M) descends 
to an identification 

5M~r\5M~r\PSL2( 



2.2.1. Automorphy conditions. Let A; > be an integer, and / : H — ^ C is a 
function on the upper half-plane satisfying the (weak) automorphy condition 

(2.1) f{^{z)) = icz + drfiz), V7=(" ^jer 

We define Ff on SM by 

Ff{z,0 :=CV(^) 

Then 

Ffh{z,0) = Ff{z,o, V7er 
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that is Ff is F-invariant so descends to a function on SM = r\S'EI, and via 
the identification F ^ F gives a F-invariant function Ff on PSL2 (M) : 



and therefore the function Ff on the group PSL2(M) transforms under the 
right action of = S'0(2)/{ib/} by the character K(a) ^ g*fc"_ 

2.2.2. Closed geodesies on M . We consider closed geodesies on M, that is 
an initial condition (-ZqjCo) £ SJl so that there is some T > and 7 G F 
with 



(the equality is in PSL2(M), that is the matrices agree up to a sign). 

Changing the initial condition (-^OjCo) to a F-equivalent one (^i,Ci) = 
<^(2^0)7o)) (5 G F (so that we get the same point in SM = F\S'H) replaces 

7 by its conjugate S"/5~^ . Thus we get a well-defined conjugacy class 7^7 
corresponding to the geodesic C. The conjugacy class is hyperbolic as its 
trace satisfies |tr7c| = 2cosh(T/2) > 2. 

2.3. A correspondence with binary quadratic forms. An binary qua- 
dratic form f{x,y) = ax^ + bxy + cy^ (also denoted by [a,6, c]) is called 
integral if a,b,c are integers, and is primitive if gcd(a,6, c) = 1. The dis- 
criminant of / is 6^ — 4ac. The modular group SL2(Z) acts on the set of 
integral binary quadratic forms by substitutions, and preserves the discrim- 
inant. 

There is a bijection between SL2(Z)-cquivalence classes of primitive bi- 
nary quadratic forms of positive (non-square) discriminant and primitive 
hyperbolic conjugacy classes in PSL2(Z) defined as follows: Given a primi- 
tive hyperbolic element 



Ffijg) = Ff{g), V7 G F 
Moreover, Ff has K-typc k since from the definition we find 

Ff{z,e^-0 = {e'"Qff{z) = e'^"F/(z,C) 




and hence that 



(2.2) 





the corresponding form is 



(2.3) 



B{l) = 



sign (a -|- d) 



[6, d — a, — c] 



gcd(6, d — a, — c) 
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which is primitive by definition, and has discriminant 
Moreover 

S(-7) = S(7), 5(7-1) = -5(7) 

Given a primitive integral binary quadratic form / = [a, b, c] of positive 
non-square discriminant d = b^—4ac, let {to, uq) be the fundamental solution 
of the Pell equation — du^ = 4 with to > 0, «o > (which exists since we 
assume d > is not a perfect square). Define the matrix 



7(/) := 



-cuo J 



which is hyperbolic of trace to = V du^ + 4 > 2, and is primitive. Then 
B{'y{f)) = f and gives a bijcction between primitive hyperbolic conjugacy 
classes in PSL2(Z) and equivalence classes of primitive binary quadratic 
forms of non-square positive discriminant. 

2.4. Periods. Consider a (primitive, oriented) closed geodesic on M; it is 
determined by a primitive hyperbolic conjugacy class 7 G F, Let C be the lift 
of the closed geodesic to to the unit tangent bundle SM. For any function 
F on SM, we define the period of F along C by choosing a point on the 
lifted geodesic (2:0, Co) (that is an initial condition) and setting 



[ F := [ Fo^*{zo,Co)dt 
Jc Jo 



where T > is the length of the geodesic, that is the first time that 
*^(^o,Co) =7(^o,Co)■ 
2.4.1. An alternative expression for the period. For a hyperbolic matrix 7 = 
^ ^ , define a binary quadratic form (not necessarily primitive) 

Q^{z) = cz^ + {d - a)z - b = z) {z - 7(2;)) 

Note that Q--y = —Q-y- 

The two zeros 'w± of are the the fixed points of 7, which are the 
intersection with real axis of the semi-circle in the upper half-plane which 
determines the closed geodesic. By (2.2), the fixed points w± of 7 on the 
boundary are 50 (0) and (70(00): Indeed, ^{w) = w iS 

g-T/2 1 % (.w) = go {w) 

that is iff e^gQ^{w) = gQ^{w), and since T ^ this forces gQ^{w) = 0, 00. 
Thus we find 

Q^{z) = C{z - gomiz - go{oo)) 



14 
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Let / : EI ^ C satisfy the automorphy condition (2.1) of weight 2k for T, 
and set F = Ff, that is F{z, Q = C,^f{z) which is a T-invariant function on 
5EI, that is a function on SM, which transforms under SO{2) with K-type 
k. Let 

rkif,l)= / fiz)Q^{z)''-Uz 
J zo 

where Zo lies on the semi-circle between the fixed points of 7 and the contour 
of integration^ is along the geodesic arc linking zq and jzq. 
Let 

:= tr(7)2 - 4 = disc(Q^) 

be the discriminant of the quadratic form Q^. Then rfc(/, 7) is simply related 
to the period of / on the geodesic defined by 7 [13, proposition 4]: 

(2.4) rfe(/,7) = (-sign(tr(7))/D;)'"' / F 

Therefore in terms of the corresponding binary quadratic form B{x, y) = 
5(7) (x,2/) (2.3), we get 

(2.5) / F = / f{z)B{l, -zf-^dz =: J(S) 

Jc (discS)~ Jzt, 

Note that the RHS above makes sense also for non-primitive forms, and is 
dilation invariant: J{tB) = J{B). 

3. Symmetry considerations 

3.1. Background on the representation theory of SL2(M). Let vr be 
an irreducible infinite dimensional unitary representation of SL2(M) on a 
Hilbert space H which factors through G = PSL2(M). Let K = SO{2), 
and let H^^^ be the space of K-&mte vectors in H, that is vectors whose 
translates by K span a finite dimensional subspace. Then 7i^^^ is dense in 
H and consists of smooth vectors, and the Lie algebra SI2 acts on H^^^ by 
dir, the differential of the action of G. 

According to Bargmann's classification of such tt's (we follow the exposi- 
tion in Lang [22] ) , there are orthogonal one-dimensional subspaces Tin , with 
n even, which are ivT-invariant and together span T^W. To be more precise, 
we consider two cases: 

i) That there is no highest or lowest K-type, this being the spherical, or 
Maass case: 

(3.1) 7^(-^)= Hn 

n even 



"^If / is holomorphic, the integral is independent of the contour 
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with Tin one dimensional for n even, say Tin = Ccpn, and the (pn satisfy 

dTr{W)(l)n = in(t>n 
(3.2) dTr{E-)4>n = (s + 1 - n)<Pn~2 

d-K{E^)(t)n = (s + 1 + n)4)n+2 

where 

(3-3) -^(i_"0. ^^(?;), -^(_», J) 

are the standard basis of the Lie algebra s[2(M), 

(3.4) E^ = H±iV 

are in the complexified Lie algebra 512(C), and dTr{E^) are the weight rais- 
ing/lowering operators. Here s G C is a parameter which, since we assume 
that TT is unitary, lies on the imaginary axis iM or in the interval (—1,1). Note 
that since wc arc assuming the representation factors through G = PSL2 (M) , 
only even weights appear. 

ii) H has a lowest or highest if -type. In the first case there is an even 
positive integer mo > so that 

(3.5) 7^W= Hm 

m=rao 

m even 

with Tijn one-dimensional, say Hm = C^^ and the satisfy (3.2) with 
s = mo — 1. In particular, 4>rno is annihilated by the lowering operator: 

(3.6) dTr{E-)4>m, = . 

These tt's correspond to holomorphic forms of even weight. 

In the case there is a highest if -type, there is a negative even integer 
mo < so that 

mo 

(3.7) 7^W= Hm. 

m=— oc 
m even 

Again Hm = C(pm for m < mo even, so that (pm satisfy (3.2) with s = 
—mo — 1 and the highest weight vector (pm^ is annihilated by the raising 
operator: 

(3.8) dTr{E+)ct)mo = • 

In case (i) we denote by the if -invariant (spherical) vector We 
normalize it so that {(j)o, (j)o) = 1 and then it is unique up to multiplication 
by a complex scalar of unit modulus. In case (ii) we denote by (p^^ the 
similarly normalized lowest/highest weight vector (pmo- We will call these 
(pT^-^s "minimal vectors" of the representation. 
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3.2. Linear forms. We consider linear forms r] on W(^) which are invariant 
by the "geodesic flow" and "time reversal symmetry" , that is 

• rj is annihilated by dTT{H), where H = ) G s[2 is the infini- 



^0 -1^ 

tesimal generator of the group of diagonal matrices A 

(3.9) ri{dTT{H)v) =0, G H^^^ 

In this case we say that rj is A-invariant^. 
1 
-1 



rj is fixed by 7r( 



(3.10) Vi7r(^\ ^^v)=n{v), Mven^^'K 

In this case we say that rj is invariant under time-reversal symmetry. 

Proposition 3.1. Let vr he an irreducible infinite dimensional unitary rep- 
resentation of SL2(M) on a Hilbert space H which factors through G = 
PSL2(M). Then the space of linear forms rj on invariant under A and 

w is at most one- dimensional, and any such form is completely determined 
by its action on a "minimal" vector cp-^. In the case (ii) of discrete series, the 
space of A-invariant forms is one-dimensional, when m = 2 mod 4 none 
of them is w-invariant, and if ni = mod 4 then any A-invariant form is 
automatically w-invariant. In the spherical case the space of linear forms 
invariant under A and w is one dimensional. 

This is shown by giving an explicit formula for ri{(l)n) in term of rj{(pT^). 
Since the cases (i) and (ii) have slightly different features we deal with them 
separately. 

In case (i), = is the spherical vector. We are assuming that rj is 
invariant under time-reversal symmetry, that is (3.10) holds. Since 

1^ 
-1 0. 



TT 



if n = 2 mod 4, due to (3.2) and oj ~ ^^P(f follows from 

(3.10) that 

(3.11) rj{(J)n) = 0, ifn = 2 mod 4 . 
Now 2H = + E~ and from (3.2) we have 

2dTr{H)cPn = {d7r{E+) + d7r{E-)) 

= (s + 1 - n)(pn-2 + (s + 1 + n)(j)n+2 ■ 

Hence 

(3.12) r]{2dTr{H)(l)n) = (s + 1 - ra)??(</'„-2) + (s + 1 + n)r/(0„+2) . 



^This choice of terminology is imprecise since n{A} need not preserve the space of 
if -finite vectors on which r] is a^priori defined 
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But the LHS of (3.12) is zero since we are assuming (3.9). Hence for n even 
and in particular n = 2 mod 4 we have 

(n - S - l)?7(0n_2) = {n + s+ l)ri{(l)n+2) ■ 
It follows that for m > 4, m = mod 4 that 

, ,, , (l — s)(5 — s) (m — 3 — s) ,, , 

(3.13) ,(0„) = M<p-n,) = ;„-i+> <^°) • 

This together with (3.11) determines r/ on T[^^^ explicitly in terms of rj^cjiQ). 

Conversely, (3.11) and (3.13) with r]{(t>o) = 1 define a unique A- and w- 
invariant linear form on TC, which we denote by (,n,4>T,- So in this case the 
space of such linear forms is one-dimensional and any such form rj satisfies 

We turn to case (ii) and show that the space of ^-invariant linear forms 
on is one-dimensional. Consider say the lowest weight case: Take the 

lowest weight vector ^mo, "t-q > even. From (3.2) and (3.6) we have 

2d'K{H)(pmo = {dTr{E+) + dTT{E')) (pmo = 2mo^!>mo+2 • 
Hence assuming rj is ^-invariant we get that 

Furthermore for m > mo even and by (3.2) we have 

(m - mo)r]{4>m-2) = {m + 777,0)77(^^+2) • 

Hence 

(3.14) vi4^m) = 0) for m > mo, m = mo -|- 2 mod 4 
and 

= („„+i)'(I„+3)'.'.''.(„!!+i-i) '''^-°' 

for = mod 4, > 4. 

Thus the space of A-invariant linear forms on H^^^ is one-dimensional. 
It is spanned by where (■Hn,4>n{(t^mo) = 1 and is defined by (3.14) and 

(3.15) . Again any A-invariant linear form rj on H^^^ satisfies 

(3-16) ri = rj{(f)^)^n^^^^ . 

The case of highest weight vectors and ^-invariant forms is the same. 

If we now impose the further condition that rj be w-invariant in these case 
(ii) representations, then invariance under time-reversal symmetry gives as 
in (3.11) that 

vi'Pm) = 0, for m = 2 mod 4 . 
This coupled with (3.14) means that if rriQ = 2 mod 4 then 77 = 0. That is 
of mo = 2 mod 4 then there is no non-zero linear form invariant under A 
and w. 
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If mo = mod 4 then from our discussion, every A-invariant hnear form 
is automatically w-invariant and in this case such hnear forms satisfy (3.16). 

3.3. Orientation reversal symmetry. We now examine the role of an 
additional possible symmetry, "orientation reversal" r. It need not act on 
irreducible representations of PSL2(M). What we do is given an irreducible 
unitary representation vr on a Hilbert space TC, we consider Hilbcrt spaces 
U which in the spherical case is the original representation H, and in the 
case of the discrete series TiJ^ where there is a lowest weight vector of weight 
m > 0, we define 

to be the direct sum of the irreducible representations with lowest weight m 
and that with highest weight — m. We write U^^^ for the dense subspace of 
K -Unite vectors in U. 

An orientation-reversing symmetry of is a unitary map r of which is 
an involution, that is 

(3.17) = J 
satisfying 

(3.18) r7r(PF) = -7r(PF)r 
and 

(3.19) r7r(£'+) = TT{E-)r . 
As a consequence of (3.19) and (3.17) we have 

(3.20) rTT{E-) = TT{E+)r . 

Moreover r commutes with the ^-action whose infinitesimal generator is 
H = ^{E'^ + E~) by (3.19), (3.20), and with time-reversal symmetry, that 
is with^ exp(f 7r(W^)) by virtue of (3.18). 

As our basic example we consider the orientation reversal involution on 
the function space L^{r\G) given by 

rf{x):=f{dx6-'), 5 = . 

The relations (3.18), (3.19) hold since for the Lie algebra elements H, V and 
W of (3.3) we have 

5W = -WS, SH = H6, 6V = -V6 . 



^The first tt is the constant 3.1415 . . . ! 
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3.4. Action of r on weight vectors. We first note that due to the com- 
mutation relation (3.18), r must reverse weights, that is 

^(t>n = Cn4>-n ■ 

with |c„| = 1 since r is unitary, and c„c_„ = 1 since = /. In particular, 
in the spherical case when there is a vector (f)Q of weight 0, we must have 

(3.21) r<t>o = e4>o, e = ±1 . 

We say the spherical representation U is even if if the sign is +, and odd if 

the sign is — . 

In the case of the discrete series representations U = 'H"^™©W~'", m > 0, 
we choose a lowest weight vector 4>m G T-C^"^ of unit length, which we call 
the minimal (or generating) vector, that r(f)m is a unit vector of weight —m, 
and normalize a choice of highest weight vector of unit length by taking 

(3.22) (f>_rn ■■= r^m ■ 

We claim that the choice of minimal vector (pQ in the spherical case and 
in the discrete series case uniquely determine r. 

Indeed, starting with 0^, the lowest weight vector for m > 0, we get from 
(3.2) for A; > 

^ k-l 

0m+2fe = -r-. 7T7r(^+) Vm, c(s; m,k) = Y[{s + m + l + 2j) 

C{^S, Ul, K) 

and for the highest weight vector (j)-m = r^m 



c(s; m, k) 



Therefore using (3.19) 



c(s; m, k) 



c(s; m. A;) 



-m—2k 



c(s; m, k) 
and likewise 

r0-m-2fe = (t>m+2k 

That is for the discrete series r exactly interchanges (pn and (f>—n' 

(3.23) r(j)n = ^-m \n\>m, n = m mod 2 

In the case of the spherical representations, the same analysis shows that 

(3.24) r(f)n = e(t)-n, n G 2Z 
where e = ±1 is determined by (3.21). 
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3.5. r-invariant functionals. If is a linear functional on U, invariant 
under the action of A and time-reversal symmetry w, then the functional 

rf : V ^ rj{rv) is also invariant under A and w since r commutes with A 
and with w. Wc wish to determine when r]^ = rj. 

Proposition 3.2. The space of linear functionals on U^^^ which are in- 
variant under A, w and r is at most one dimensional. In the spherical 
case there are no such functionals for odd representations, and the space is 
one- dimensional in the even case, every functional invariant under A and 
w being automatically r-invariant. For the discrete series there are no such 
functionals for weight m = 2 mod A, and for weight m = mod 4 the space 
of A-invariant functionals is two-dimensional, each is automatically invari- 
ant under w and the suhspace of r-invariant functionals is one-dimensional. 

Proof. We start with the spherical case. There is a one-dimensional space of 
functionals invariant under A and and we take the unique one satisfying 

riiM = 1 • 

Hence rf , being itself invariant under A and w, must be a multiple of rj, and 
because = 7 we have 

rf = ±rj . 

We claim the sign is determined by the sign in (3.21), that is if r^o = e^o 
then 

rj'^ = erj . 

It suffices to check this on the spherical vector (j)Q, that is to show r]'^i4'o) = e. 
Indeed, we have 

'?'"(0o) = ?7(r0o) = vi^iM = ^r]{(t)o) = e 

as required. Thus in the odd case rj^ = —r) and there are no A- and w- 
invariant functionals which arc r-invariant, and in the even case rf = r] and 
every A- and w^-invariant functional is automatically r-invariant. 

In the discrete series case, there are no functionals invariant under A and 
w if m = 2 mod 4, hence we only consider the case m = mod 4. In that 
case there are unique A-invariant functionals 77+ on 7^+"^ and rj- on 'H~'^ 
satisfying 

??+(0m) = 1, V-i(P-m) = 1 

and these are automatically invariant under time reversal symmetry. Hence 
the space of A-invariant functionals onU = 'H'^'^ ®'H~'^ is two dimensional, 
consisting of linear combinations 

r] = c+??+ © c„7?_ 

and these are automatically invariant Tinder time reversal symmetry. They 
are uniquely determined by their action on the lowest and highest weight 
vectors and ^-m = ^(l>m- 

c± = rji<P±m) 
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Since t]^ is also yl- invariant, we have 

r?"" = c'_,_r/+ + c'_r7_ 

Now T]'^ = 7] ii and only if = c+ and c'_ = c_. We claim that this 
happens if and only if c+ = c_ , which will show that the space of A-invariant 
functionals which are r-invariant is exactly one-dimensional in this case. 
Indeed, we have 

c+ = r]''{(pm) = r]{r4>m) = vi't'-m) = C- 
c'_ = rf{(t)-m) = V{^4>-m) = V{<Pm) = C+ 

and so c'_|_ = c± if and only if c+ = c_ as claimed. □ 

3.6. Bilinear forms. Wc apply the above uniqueness of linear forms to bi- 
invariant sesqui-linear forms onU xU. Let T{v,v') be such a form, that is 
linear in v, conjugate- linear in v' and invariant under A, w and r in each 
variable separately. For instance, we can take 

J 

T{v,v') = Y,rij{v)r,'j{v') 

where rjj, are invariant linear forms. From the prior discussion, 

T{v,v) = 

if TT is of type (ii) with mo = 2 mod 4. Otherwise T is completely deter- 
mined by value T{^-m 4>-k) at the minimal vector ^^r- In feet T is the product 
of linear forms 



T{v,v')=T{cl)^,4>^)^u,^^{v)^u,^^{v') . 

where S,u,<f>T, is the unique invariant linear form taking value 1 at the minimal 
vector 

3.7. Application to the classical and arithmetic variances. We apply 
these remarks to the measures and to the classical variance V. We 
consider the discrete decomposition of the regular representation of G = 
PSL2(M) on L^gp(r\G). For an irreducible sub-representation, form the 
space as above. 

3.7.1. The arithmetic measure ji^ is a linear form on 14-,^ invariant under A, 
w and r. Hence ^d{F) = if vr is a discrete series with weight mo = 2 
mod 4 and otherwise 

(3.25) ^id = fJ'd{(l)7v)£,u^,4>^ ■ 

Hence if Fi and F2 are in Z^tti and W,r2 the sesqui-linear jid sums take the 
form 

E MFl) li-djFi) 7^ /id(07riW(07r2) 

d<Y Vd ^ d 
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This gives a universal reduction for computation of the variance of jid to the 
cases Fi = (f>Tr^ , F2 = (f>TT2 ■ 

3.7.2. The classical variance V is by its definition diagonalized by the irre- 
ducibles in the decomposition of L^{r\G). We define projections onto the 
set of w-invariant functions 

jpeven ._ i(F + F^) 

and onto the set of functions invariant under both w and r, 
Set 

We wish to completely determine V^^^ and V^'"^'^. 

For an irreducible vr, V^^^ vanishes on Ut^ if tt is discrete series of weight 
rriQ = 2 mod 4 and otherwise is given by 



It remains to determine ^'^^'^((^Tr, '/'tt)- 

Lemma 3.3. i) For tt spherical with parameter s = ir, 

(3.26) y.>»(^„,^) = _I_LV-JL(^,^) 

a) For TT discrete series of weight m = mod 4 (m > 0), 

(3.27) V'y^{(f>m,(f>m) = 22"^Beta(-,-) 
where Beta is Euler's beta function. 

Proof. In the spherical case we need to compute V^y'^{(f)Q, 00 ). For spherical 
representations, wc saw that yl-invariance and w-invariance automatically 
imply r-invariance, hence on such spherical lA, 

Moreover, since (^o is spherical, w^o = 00 and hence ^q^"^" = 0o- Thus 
F*f"^(0o,0o) = ^(00)00)) which was computed in [25] to give 

l^--(0o,</>o) = V{Mo) = 2j|r4 + 2ir)P ^^°'^°^ • 

For TT discrete series of weight mo = mod 4 (m > 0) then ^-invariance 
implies w-invariance, hence 

T reven \ 

V \UxU = V \uxU 

and so 

^'^"*(0m, 0m) = \ {Victim, <t>m) + ^(0m, r0^) + ^(r0m, 0m) + ^(r0m, r0^)) 
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By its definition, V respects the orthogonal decomposition into irreducibles; 
since (f)m G TTm and rcf)^ = ^-m € T^-m lie in distinct irreducibles, we get 
V{(f)Tn.,v(f)Tn) = = V{v(f)m-,<t)m)- Moreover we have 

V{vFi,vF2) = V{Fi,F2) 

for any Fi, F2. To see this, note first that r is induced by the measure 
preserving map x 1— > 5x5~^ of SX = r\G and hence 

{rFi ,rF2)= [ rFi {x)vF2{x)dx = {Fi , F2) 
Jt\g 

Moreover, r commutes with the geodesic flow and so 

V{rF,,rF2) = J^Jtt ("^^q^' ^_%) m,vF2)dt 

= V{FuF2) . 

Thus we find 
Let 

/(x) = {7r{x)(j)-m, (p-m), X £ G . 

Then applying the raising operator E'^ via the regular representation gives 
an operator >C^+ satisfying 

(3.28) jCE+f{x) = {TT{x){d'K{E+)(l>_^),(l)_m) = 

Also by the unitarity of tt, 

f{k{ei)xk{92)) = {n{k{Bi)xk{B2))<t>-mA-m) 

\ 

Using the coordinates k{6i) I ^ -r/2 1 ^(^2) on G and the formula for 
in these coordinates, we deduce from (3.28) and (3.29) that 

where g satisfies the ODE 

,„ do coshr m 

3.30 2/ = -^-r-m^ + —-g 

dr sinh r smh r 

and since we have normalized ((^—j^, 0—^^) — 1, we have 

5(0) = 1 . 
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We integrate (3.30) and find that 



g{r) = (cosh^)" 



Hence 



(cosh-)— dr = 2-Beta(-,-) 



-oo 

where Beta is Euler's beta function. Thus we find 

1 TJl 771 



4. Half-integral weight forms 



□ 



4.1. Basic properties. Let F be a discrete subgroup of SL2(M) of finite co- 
volume. Given a character x : T — >^ S-*^, an automorphic function of weight 
k and character % for F is a function / : H — C satisfying 

with suitable growth conditions at the cusps of F. It is cuspidal if it vanishes 
at the cusps. 

The Laplacian of weight k is defined as 

The Laplacian A;;, maps forms of weight k to themselves, and maps cusp 

forms to themselves. A Maass cusp form of weight /c is a cuspidal automor- 
phic function of weight k (for some character x) which is an eigenfunction 
of Afc. 

Let Wi^^fj, be the standard Whittaker function, normalized so that at in- 
finity 

(4.1) w^«,/.(y) ~ y'^e"^/', 2/^00 

The asymptotic behaviour of W^^^iy) near y = is 

F(|-/x-k) T{^+ij.-k) 

for |Lt 7^ 0, while 

(4.3) VF„,o<y'/'logy, y^O 

The functions 

giz,s) :=W^|_^,_i(47ry)e(±x) 
are eigenfunctions of with eigenvalue A = s(l — s). 
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A Maass cusp form F of weight k and eigenvalue A = — s) has Fourier 
expansion 

Fiz) = ^p(n)/f "(")(|n|z,s) = T P(n) W-.i,n(.). ^ , (47r|n|j/)e(nx) 

The Petersson inner product is defined for a pair of (cuspidal) functions 
of the same weight A; and character as 



-dxdy 

\H " ' ■ " ■ y'^ 



if, 9) = [ f{z)g{z) 

Jr\m 



4.2. Maass operators. For any real A;, define the raising operator 

k f.d d\ k , ,d 

and the lowering operator 

k / d d\ k , , d 

The raising operator takes Maass forms of weight k to forms of weight 
k + 2 and the lowering operator takes Maass forms of weight k to forms 
of weight k — 2. 

Then 

The effect of the Maass operators on Petersson inner products is given as: 
If /, g have weight A; and character %, then 

(Kfc/,M) = rA(5)-A(-M {f,g) 



and 



(Afc/,Afc<7) = (A(5)-A(^))(/,5) 



The action of the Maass operators on the eigenfunctions f^{z,s) is 

(4.4) Kkf^{z,s) = -f+^,{z,s), Kkf^{z,s) = {s+'^){l-s+'^)f^^,{z,s) 
and 

(4.5) A,f+(z,s) = -{s-'^){l-s-'^)f^_,{z,s), Akf-iz,s) = f^_,{z,s) 
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4.3. Maass operators and Fourier expansions. We want to see the 
Fourier expansion of a "raised" Maass form in terms of its original. 

So start with a Maass form F of weight 1/2 and eigenvalue A = 1/4 + 
with Fourier expansion 

^i^) = T^P(")^ Bign(n) ..^(47r|ra|ii/)e(raa;) 

Applying the Maass raising operator S^t a form of weight 

5/2 whose Fourier expansion is obtained by (4.4) as 

oo 

(4.6) Ki/2F(z) = -p(n)W^5/4,ir(47rny)e(nx) 

°° / 3 \ 
+ ( (4)' + ) P(-n)^-5/4,ir(47rny)e(-nx) 

n=l ^ ^ 



n=l 



Applying the lowering operator we get a form /^i/2F of weight —3/2 
with Fourier expansion is obtained by (4.5) as 

(4.7) A,/2F{z) = f; - U^)^ + r^") p(n)VF_3/4,^.(47rny)e(nx) 

n=l ^ ^ 

00 

+ Yl P(-"')^3/4,ir(47rn?/)e( 



n=l 



5. Rankin-Selberg theory 

5.1. Classical Rankin-Selberg theory. We recall classical Rankin-Selberg 
theory as applied to a holomorphic form F of weight fc + 1/2 with Fourier 
expansion 

d>l 

Let E{z, s) be the standard Eisenstein series for ro(4): 



where 



T6roo\ro(4) 



The series is absolutely convergent for 3FJ(s) > 1, with an analytic contin- 
uation to 3f?(s) > 1/2 except for a simple pole at s = 1, where the residue 
is 

^^^-^^^"'^) = vol(ro(4)\M)=i 



VARIANCE OF ARITHMETIC MEASURES ON THE MODULAR SURFACE 27 

One starts with the integral 

/(.)= / \F{z)\'E{z,s)y'^+'^^ 
Jro(4)\H y 

which is analytic in 3?(s) > 1/2 except for a simple pole at s = 1 with residue 

By the "unfolding trick" , we have 

/(.) = (4.)-(^+'=-|)r(. + k-hf: l^pn- 

n=l n 2 

and hence the Dirichlet series 



CF{n) |2^_, 



fc-l/2 

n=l n 2 

has a simple pole at s = 1 with residue 

(4^)^+^ {F,F) 

Consequently we find 

1 v-,CF(n),2 (47r)^+^ {F,F) 



(5.1) lim ^ y 



Similar considerations show that if we take forms F of weight A; + 1/2 and 
G of weight ^ + 1/2 {k and £ possibly different) which are orthogonal then 
we have 



(5-2) 1 V ^ 



CF{n) coin) 



n<N n 2 n 2 







and that if F is a Maass form f weight 1/2 for ro(4) with Fourier expansion 
F{x + iy) = ^ p(n)VFsign(„)/4,jr(47r|n|y)e(na;) 

and G is a holomorphic form of weight A; + 1/2 then 

(5.3) J^^^E^^H^ = 

n<N n 2 

These arguments will also give the asymptotics of the sum of squares 

^ |47rnp(n)|2 
-Ar<n<Ar 
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of Fourier coefficients with both positive and negative indices. However 
for our apphcation we need to be able to separately sum only coefficients 
indexed by positive integers, that is we require the asymptotics of the series 

N 

^|47rnp(n)|2 

n=l 

To do so, we make use of the arguments in the paper by Matthes [28] , which 
we adapt for our case, see also [27, 29]. 

5.2. One-sided Rankin-Selberg theory. Let F{z) and F'{z) be Maass 
cusp forms of weight 1/2 for ro(4), and Laplace eigenvalues A = ^ + r^, 
A' = I + (r')^, with Fourier expansion 

(n)/4,ir i^'^ I ^ 1 2/ ) ^ ) 

F'{x + iy) = ^p{n)W^i^(^ri)/i,ir'{^An\y)e{nx) 
We define two Dirichlet series 



n=l 



oo 



n=l 



Proposition 5.1. Both L±{s) have analytic continuation to K(s) > ^, 
except for a simple pole at s = 1 if F and F' are not orthogonal. 

We next compute the residue at s = 1 when F' = F: 

Proposition 5.2. The residue at s = 1 of L^{s,F x F) is 

(5.4) R-^ := ReSs=,L^{s,F x F) = ,1^/,^^ 

and the residue of L-{s,F x F) is 

|r(|+^r)|2 {F,F) 



R~ = 



r(i + 2zr)|2 TT 



The arguments and ideas needed to establish this have been essentially 
provided in [28]. 

As a consequence, we deduce by a standard Tauberian argument that 
Corollary 5.3. Let F be as above. Then 
(5.5) 47rn|p(n)|2 ~ R+N, N ^ oo 

l<n<N 
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while if F and F' are orthogonal then 

(5.6) A'Knp{n)p'{n) = o{N) 

l<n<N 

Proof. Applying the Wiener- Ikehara Tauberian theorem (see [23] for exam- 
ple) to the Dirichlet series 

E°° 47rn|p(n)p 47rn|/9'(n)p ^ 47:n\p{n) + p'{n)\'^ 

n=l n=l n=l 

respectively, we infer by proposition 5.1 that, as N ^ oo, 

^ 47rn\pin)\^ ^ R+N, 

l<n<N 

4^n|p'(n)|2~i?+Ar, 

l<n<N 

and 

Y 47rn|p(n) + p'(n)|2 ~ {R+ + R+)N. 

l<n<N 

Thus, 

Y 4™ ^{p{n)p'{n)) = o{N). 

l<n<N 

Similarly applying the Wiener- Ikehara theorem to the Dirichlet series 

A-!Tn\p{n) + ip'(n)p 



E 

we obtain 



n=l 



J2 47rn|/9(n) + ip'{n)\'^ ~ {R+ + i?+ )^, 

l<n<iV 

and consequently, 

Y 4vrn 9(p(n)p'(n)) = o(iV). 

l<<iV 

This proves 

i7rnp(n)p'{n) = o{N). 

l<n<N 

□ 
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5.3. A Mellin transform. Fix fi = ir and ly = ir' and define 

PCX) 

Mkis) := / WkMWkAy)y'~^dy 

Jo 

In view of the asymptotics (4.1), (4.2), (4.3), the integral is absolutely con- 
vergent for ??(s) > and hence A^fe(s) is analytic in that region. The 
asymptotic behaviour of 7Wfc(,s) is given by 

Lemma 5.4. Assume that \k\ < 1/2. Then as \s\ oo, 

A4.(.) = ^^^|i^(i + o,,(H-/^)) 

Proof. This is a direct generalization of Lemma 4. 1 in [28] , which deals with 
the case r' = r. As in [28], we use the integral representation 

If T{l-v-ir)r{^-v + ir)T{v-k) 

^TTZ 7l 1 (2 + «r — A;)r(2 — ir — k) 

where the the path of integration runs from —ioo to zoo and is chosen so 
that all poles of T{v — k) are to the left, and all poles of r(^ — v ±ir) are 
to the right of L; this is possible since we assume that \k\ < 1/2. Inserting 
this into the formula for Mk{s) gives 

Mk{s) = 

Jo 27rz 7i +ir - k)r{^ - ir - k) 

Then one uses the formula 

to find 

Mk{s) = ^ I T{-l + v + s + ir')r{-l+v + s-ir') 
2m Ji^ 2 2 



r(^ - t; - ir)r(i -v + ir)r{v - k) 



-dv. 



r{v + s- k)r{^ +ir- k)r{l -ir-k) 

One then shifts the contour of integration to the line 3f?('u) = k — 1/2, picking 
up a single residue at v = k, and estimates the remaining integral as in [28] 
giving 

T{s-l + k-ir')T{s-^ + k + ir') ^ ^ {^^ 3 



Mk{s) = ^ + O^,, [V{s - 2 + 2fc) 

The conclusion of the Lemma now follows from Stirling's formula. □ 
Let 

~ M-k{s + l) M-k{s) 



(5.7) M{s) 



Mk{s + l) -Mkis) 
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Lemma 5.5. M(s) is analytic and nonsingular for 5}?(s) > 0. 

Proof. Holomorphy in 3fi(s) > follows from that of M.k{s). As in [28], one 
shows that there is a recurrence relation 

(5.8) (s + l)Mk{s + 2) - 2k{2s + l)Mk{s + 1) 

By the recurrence relation (5.8), we infer that 
(. + 1) det Mis + !) = (. det M(.)) " + ^^^^ " " '^^'^ , 

and 

r(s + n + i^)T(s + fj- — '^)r(s — ^ + i^)^is — — y) 



det M(s) 



r(s + n + /X + i')T{s + n + n — v)T{s + n — /x + v)T{s + n — jx — v) 
(s + n)r2(s + n) 



■ det M{s + n) 



s r2(s) 

By using Lemma 5.4 with Stirling's formula, we deduce that^ 

(s + n)r2(s + n)detM(s + n) 
n^oo r(s + n + 11 + z^)r(s + n + — z^)r(s + n — + i')T{s + n — ji — u) 

= -2 . 

Therefore we conclude that 

r(s + IJL + i')T{s + 11- i')T{s - IJL + z^)r(s - 11- u) 



det M(s) = -2- 



sr2(s) 

and thus det M(s) / for 3fi(s) > 0. □ 

5.4. Proof of Proposition 5.1. We define the Eisenstein series of weight 
—2 and 2 by 



7er=o\r(,(4) 



These series are absolutely convergent for 3fi(s) > 1, with an analytic con- 
tinuation (no poles) to 3ft(s) > 1/2. 



'Note the limit is —2 instead of 1 as in Matthes' paper [28]. 
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Using the Maass raising operator Ki/2) we get a form of weight 5/2 

with Fourier expansion 

oo 

^i/2F{z) = -p(^)^5/4,ir(47rny)e(nx) 

n=l 

n=l ^ ' 

and using the lowering operator A1/2 we get a form A1/2F of weight —3/2 
with Fourier expansion 

Ai/2^(^) = X] - ((\)^ + P(n)VF-3/4,ir(47rn2/)e(nx) 

n=l ^ ' 

00 

+ P(-"')^3/4,ir(47rn?/)e(-na;) 

n=l 

Consider the Rankin-Selberg integrals 

dxdy 



A{s)= [ F{z)F'{z)E{z,s)- 



/ro(4)\H 
and 

dxdy 



y2 



2B{s)= [ Ky2F{z)F'{z)E^2{z,s) 

"'ro(4)\H 

+ / A F{z)F^)E2{z,s)^ 
Jro(4)\H y 

Since E±2{z,s) are analytic in 3ft(s) > 1/2, B{s) is analytic in 3ft(s) > 1/2; 
and since E{z, s) is analytic in > 1/2 except for simple pole at s = 1, 
A{s) is analytic for 3fJ(s) > 1/2 except for possibly a simple pole at s = 1, 
where the residue is 

(F F') (F F') 

Res i^fs) - ^ - ^ 

so that ^(s) is analytic also at s = 1 if F and F' are orthogonal. 
By the standard unfolding trick, we find that for 3ft(s) > 1, 
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and 



4:7rnp{n)p'{n) 



n=l 



{47m) 



n=l 



Setting k = 1/4 we then have 

A{s) = (47r)-^ {L+{s)Mk{s) + L-{s)M-k{s)) 

and moreover 
Lemma 5.6. 

(5.9) B{s) = {47r)-'{L+{s){r/2Mk{s) - l/2Mk{s + 1)) 

+ L^{s){r/2M-kis) + l/2M-k{s + l))) 

Proof. Let k > he any half integer. In our case k = 1/2. We have, by 
unfolding the integral, that 

/ KkFiz)F^E^2iz, s)^ = r C KkF{z)F^)f-''dxdy, 
JTo\n y Jo Jo 



and 



AkF{z)F'{z)E2{z, s)^= r f\kF{z)F'{z)y'-^dxdy. 

[ y Jo Jo 

Now 

KkWk^g^^^)^.^i4TT\n\y)e{nx) 



ro\H 



{k/2 - 2sign(n)7r|n|y)VF|^^„(„)_ i,.(47r|n|y) 
+ 4Tr\n\yW,^^^^ JMn\y)] e{nx), 



and 

^kW.^^^.^{4Tr\n\y)e{nx 



□ 



{k/2 - 2sign(n)7r|n|y)W|^^„(„)_ .^(47r|n|y) 
- 4Tr\n\yW',^^^^^^^ ,,(47r|n|y)] e{nx). 

Hence Lemma 5.5 follows. 

Consequently, 

A{s) = (47r)-*- {L+{s)Mk{s) + L_{s)M-k{s)) 
rA{s) - 2B{s) = (47r)-*- {L+{s)Mk{s + 1) - L_{s)M.k{s + 1)) 
Solving for L^{.s), we obtain 



L+{s) 
L-{s) 



{47T) 



-M{s) 



A{s) 

rA{s) - 2B{s) 



det M(s)' 

where M{s) is given by (5.7). Therefore Proposition 5.1 follows from Lemma 5.5. 
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5.5. Proof of Proposition 5.2. Consider the Rankin-Selberg integrals 

dxdy 



I{s)= f Fiz)F{z)E{z,s)- 



J{s)= I Ky^F{z)F{z)E_2{z,s)- 

Jro{4)\H 



/ro(4)\H 
and 

^dxdy 

y2 

Since E_2{z, s) is analytic in > 1/2, J^{s) is analytic in 3fi(s) > 1/2; and 
since E{z, s) except for simple pole at s = 1, 1{s) is analytic for 3fi(s) > 1/2 
except for a simple pole at s = 1, where the residue is 

By the standard unfolding trick, we find that for 3fi(s) > 1, 



|2 rco 

and 

(5.11) j{s) = - f; r w,f,,Uy)w,/4,iriy)y'~'dy 

n=l ^ 

Denoting by the residue at s = 1 of L±{s), we find from (5.11) that 
since J{s) is analytic at s = 1 

((^)' + r-^) w^-5/4,i.(y)w^-i/4,i.(y)y i?- 

= /"^5/4,^r(y)^^l/4,^r(y)-^+ 

jQ y 

Prom the formula [6, p. 858] that for |3fi(/x)| < 1/2, 

1 1 



r(i/2-K + //)r(i/2-A-/Lt) r(i/2-K-//)r(i/2-A + At)J ' 

it follows that 

dx 1 2r 



POO 



X sinh(2r7r) |r(l/4 + zr)|2 ' 
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^0 



and 

{It 1 2r 

W 5 . ixW 1 . ix)— = = - 

-t'*'-^^ -i'*'-^^ X sinh(2r7r) |r(7/4 + ir)|2 ' 

Hence 

r(l/4 + ir)|2 
r(3/4 + ir)|^ 

Computing the residue at s = 1 of T{s) using (5.10) gives 

We have for |sR(/x)| < 1/2 [6, p. 858] : 

Jo > z sin27r/x r(l/2 + /x-K)r(l/2-/x-K) 

where il>{x) = ^logr(a;). Therefore 

2dy TT tl;{l/4: + ir) — 'ip{l/4 — ir) 



y isinh27rr |r(l/4 + ir)p 

and 

2C?y _ vr V(3/4 + ir) - ^^(3/4 - ir) 



/"OO 

(5.15) / W_i (yy^ 

Jo ■* 



y isinh27rr |r(3/4 + ir)p 

Inserting (5.14), (5.15) and (5.12) into (5.13) gives 

TT 



(5.16) 2{F,F) = R+ 



isinh(27rr)|r(i + ir) 



(113 3 
+ ir) - - ir) + i/ji^ + ir) - V'(^ - ir) 

Taking the logarithmic derivative of the reflection formula 



r{x)T{i-x) = 



TT 



sin(7ra;) 
gives 

^/J{x) — — x) = —IT cot{iTx) 

Hence we find 

113 3 
ip{- + ir) - tpi- - ir) + V'(^ + ir) - tpi- - ir) 

TT COS tt{j + ir) ^ TT COS tt{j — ir) 



sin7r(|+ir) sin7r(| — ir) 
TTi sinh(27rr) 27ri sinh(27rr) 



sin7r(4 + ir)P cosh(27rr) 



2isinh(27rr)|r(- + 2ir)|^ 
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Therefore we get 

|r(4 +^nr 

which proves (5.4). 

6. Theta lifts and periods 

We summarize the results on theta hfts of Shintani [41] and Kohnen 
[15, 16] in the holomorphic case, and of Katok and Sarnak [14] in the Maass 
case. 

6.1. Holomorphic forms. For a holomorphic form / of weight 2k (we 
will later take only even k) for the full modular group SL2(Z), Shintani 
[41] showed that it can be lifted to a cuspidal Hecke eigenform 9{f,z) G 

5'fc+i(ro(4),Xfc) of weight k + ^ with character Xk{(^^ ^j) = i^)'' for 
To (4), that is transforming as 

Fijz) = XkilWl, zf^+^F{z), 7 G ro(4) 
where Ji{^,z) = Oii^z) /ei{z), 6i{z) = T.n=^oo<n'' z). Thus 



with = 1 if d = 1 mod 4 and = z if d = 3 mod 4. In particular 

Ji(7,z)2 = (^)(cz + d). 

We write the Fourier expansion of 6{f, z) is 

e{J,z) = Y^Cf{d)e{dz) 

d>l 

Then 



Cf{d)= y / fiw){a-bw + cwy 

.Jc(a) 



disc(g)=d'^'^('^) 

the sum over all SL2(Z)-equivalence classes of binary quadratic forms q = 
[a, b, c] of positive discriminant d = — 4ac. 
We note that from (2.5) 

f{w){a-bw + cw'^)^~^dw = d^ / fdnc(q) 

IC{q) Jc(q) 

is a simple multiple of the period of / over the closed (primitive, oriented) 
geodesic C{q) corresponding to q. Thus 

(the sum over all forms of discriminant d, not necessarily primitive) , that is 

Cf{d) Md(/) 



(6.1) 
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6.1.1. An inner product formula. Denote by 'S'^_)_i/2(^o(4)) the space of cusp 
forms of weight A; + 1/2 on ro(4), whose Fourier expansion X]n>i c{n)e{nz) 
satisfies the condition c(n) = unless {—l)^n = 0, 1 mod 4. It was proved 
by Kohnen [15] that the two spaces 5'2fc(ro(l)) and isomor- 
phic as modules over the Hecke algebra under the Shimura correspondence. 
Assume k is even, let 

CXD 

f (z) = Y, a f{n)e{nz) e S2k{r oil)) 

n=l 

and normalize the L-function of / as 

°° a/(n) 



L{sJ) = Y. 



„5+{2fe-l)/2 
n=l 



so that the functional equation is s i— >^ 1 — s. Let 

hf(^) =^Chfin)e{nz) G 'S'+_i/2(ro(4)) 

n>l 

corresponds to / as above such that {hf,hf) = 1. By the works of Kohnen 
[16, Theorem 3 and Corollary 1], we have 

(_l)fe/22fc 

Chf{m)chf{l) = —j^j^af{m) , 

and 

Hence we see that 9{f, z) 7^ if and only if L{^, f) / 0, and in this case hf 
is proportional to 6{f, z): 

(_l)fc/22fc 



Chfil)hfiz) = ^jj-j^9{f,z) 

Thus we get a formula for the inner product of the lift: If g is another 
cuspidal Hecke eigenform in S'2fe(ro(l)), then 

(6.2) {e{f,-), e{g,-)) = i^^L{^J){f,g) 

6.2. Maass forms. Given an even Maass form (p for SL2(Z), with eigen- 
value X = 1 + (2r)^, the theta lift F{z) = 9{(f), z) is a Maass form for ro(4) 
of weight 1/2, that is transforming under ro(4) as 

F{jz) = J{j,z)F{z), 7Gro(4) 

where J{j,z) = e{jz)/e{z), 9{z) = y^'^9i{z) = y^/*T.n=_^e{n^z). More- 
over F is an eigenfunction of A1/2 with eigenvalue | -|- r^. 
The Fourier expansion of F is given by [14] as 

F{u + iv) = Vp((i)I^.ig„w UTr\d\v)e{dv) 
d/0 



38 



WENZHI LUO, ZEEV RUDNICK AND PETER SARNAK 



where for d > 

the sum over all SL2 (Z)-equivalence classes of binary quadratic forms q = 
[a, h, c] of positive discriminant d = b'^—Aac and C(g) is the closed (primitive, 
oriented) geodesic corresponding to q. (For d < there is an analogous 
formula involving Heegner points). Thus 

Moreover if -0 is another Hecke-Maass eigenform, we have, in view of [14, 
formula (5.6), p. 224], the inner product formula 

(6.4) (0(0, •), 0(0,-)) = ^A(l/2, V), 

where 

w ,^ -e-.^ ^ s + 2ir \ ^ f s — 2ir\ ^ , 

A{s, 0) = TT [-^J r [-^—) ^(^' 

Note the above formula is still valid even if 0, -ip are not even. 

7. Proof of the main Theorem 

7.1. Expected value of fi^i. In section 6 we identified the measures Hdif) 
with Fourier coefficients of theta-lifts up to a normalization. Hence the 
vanishing of the mean value of iid{f)/d}^^ follows from the corresponding 
fact for Fourier coefficients of forms of half-integer weight. 

The first statement of Theorem 1.1 follows immediately from the Hecke's 
bound that for any a G R and e > 0, we have 

ain)e{an) = 0{N^/^+') , 

n<N 

where a{n) is the normalized Fourier coefficient of any holomorphic or Maass 
cusp form, and the implicit constant depends on the form and e alone. For 
the proof, see for example [9, page 111, Theorem 8.1] and [10, page 71, The- 
orem 5.3]. One makes use of the following formula ([6, page 857, 7.611]: 



fOO 



7r3/22^' sec(//7r) 



/•c 

•''0 1(4-2 + 2)^-4 2 2- 



k4 2^2 / V4 2 2) 

7.2. Proof of Theorem 1.2. We wish to compute the bilinear form 



- ^^^00 #{d:d<N} ^ dV2 

- d<N 



where f.g G L'^^gp(T\G) are smooth and i^- finite, and the sum is over 



discriminants. We take / and g to lie in the irreducible subspaces U-,^ defined 

-"cusp ( 



in (1.11), that is subspaces of {T\G) irreducible under the joint actions 
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of G, the orientation reversal symmetry r and under the Hecke algebra. We 
wish to show that for such f,g, 

(7.1) B{f,g) = Q 

if /, g lie in distinct (hence orthogonal) subspaces Uf^ Ug, and to compute 
-B(/, /). By the results of § 3, it suffices to consider "minimal", or generating 
vectors, that is to consider holomorphic forms or Maass forms. In particular 
we need to show that for such /, 

where 

(S/tt, / Maass form 
Ii/tt /holomorphic 

Since both B and vanish when / is holomorphic of weight = 2 

mod 4 or an odd Maass form, it suffices to treat the cases of holomorphic 
forms of weight divisible by 4 and of even Maass forms. To do so, we recall 
that in these cases we have identified iJ,d{f) with simple multiples of the 
Fourier coefficients of theta-lifts 9{f,z). Thus we may use Rankin-Selberg 
theory (Corollary 5.3) to recover (7.1) and (7.2) once we have made the 
correct identifications. We treat separately the case of holomorphic forms 
and Maass forms. 

7.2.1. Holomorphic forms. For a cuspidal Hecke eigenform of weight 2k, k 
even, the theta lift has weight A; + 1/2 with Fourier expansion 9{f,z) = 
^d>i(^fid)e{dz) with (6.1) 

Cfid) _ Hdif) 



|2 -s 

n 



By standard Rankin-Selberg theory (5.1) 

(4.)-(^-^'=-l)r(. + ^-l)f 1^ 

n=l n 2 

has a simple pole at s = 1 with residue 

i?+(/) = i^i^MZl) = i^(47r)-^r(A:)L(i, /)(/,/) 
by (6.2) and hence the Dirichlet series 

n=l 

has a simple pole at s = 1 with residue 

nk^) ^^^"0Fr(A: + i)^^2'^)^^'^^- 
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Note that by (3.27) 

^-(/J)42-lM(/,/) = .^=I|L_(;,;) 

and therefore the residue at s = 1 of D{s) is 

-y^f-(/,/)L(i /). 

TT 2 

Consequently we find 

n<N 

7.2.2. Maass forms. Let (p be an even Maass Hecke eigcnform with Laplace 
eigenvalue 1/4 + (2r)^ and F = 6{(p, •) its theta-lift, with Fourier expansion 

F{u + iv) = V p(d)VF .ign(d) . {An\d\v)e{dv) 
Recall (6.3) that for d > we have 

Thus by Corollary 5.3 we have 

2 



l<n<W 

where R'^{(j)) is given by (5.4). 
Inserting (6.4) into (5.4) gives 

We note that by (3.26) 

\r(l + i^)|4 



^2 

and hence 



27r|r(i + 2ir)|2 



(7.4) R+ = l^v^yrn^^^^^L{l,cl>) 

Therefore we get 

7.2.3. Orthogonality. Finally, the fact that B{f,g) = if the subspaces Uf, 
Ug are distinct follows from standard Rankin-Selberg theory when at least 
one of f oi g is holomorphic, see (5.2) and (5.3), while the case of both /, 
g being Maass forms follows from Corollary 5.3. 
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